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ABSTRACT 


The analytical considerations which support the computation of 
(i) two components of the velocity vector from an x-array and (ii) the 
transverse vorticity from the x-array and an adjacent parallel wire pair 
are presented herein. The electronic circuit which will execute these 
computations at a 50 khz rate is also described. The factors limiting 
the accuracy of the measurements are identified and quantitative estimates 
are given for a typical x-probe. An extensive analysis of the factors 
which effect the output and which are unknown {or unknowable) during 
the measurement is presented. Quantitative estimates of these effects 
are developed in the form of an uncertainty analysis. Numerical values, 
calculated using the analytical structure of the response equations, are 
tabulated; other estimates which require special experiments are described. 



1 . INTRODUCTION 


The vorticity, which is a measure of the fluid particles rotation 
rate about its centroid, is a continuous function in space defined by the 
operation 

^ ^ V (1) 

The operation of forming the spatial derivative implies that contains 
less information than the original velocity field, V . For example, the 
vorticity is the antisymmetric part of the second order tensor which 
represents the spatial derivative of the velocity fields, specifically 
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where A is the rate of strain tensor and B is the rate of rotation tensor* 

The loss of information is, however, often compensated by an enhanced 
understanding of the flow field when the vorticity, and its time dependent 
behavior, are used for diagnostic purposes* Some examples and reference 
to several general examinations of vorticity considerations are provided 
below. 

Lighthill [1963] presents an excellent summary of vorticity funda- 
mentals and their relationship to the factors of interest in fluid mechanical 
descriptions. Vorticity in turbulent flows is given considerable attention 
by Tennekes and Liumley [1972]. The highly instructive photographs of 
large vortical structures in a shear layer by Roshko and co-workers (see, 
e* g* , Brown and Roshko [1 974]) is but one reference which suggests that 
large scale vortical motions are of critical' importance in the description 
of turbulent shear flows. Direct measurements of the vorticity which 
constitutes the elementary ingredient of the large structure is clearly 
desirable, especially as it might provide some direct experimental evidence 
in relationship to the recently advanced theoretical (and provocative) 
arguments of Moore and Saffman [1975]. Turbulent shear flows, which 
are not bounded by a solid surface, are bounded by a thin region termed 
the viscous super layer. The rationale for the existence of this layer and 
theoretical considerations for many of its properties are developed in a 
comprehensive study by Cor r sin and Kistler [1955]* The final stages of 
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entrainment of ambient flow into the shear flow is known to occur at this 
boundary; the relative importance and a description of the behavior of 
engulfment or the convective action of the pressure field associated with 
the corrugations of this interface would clearly be aided by an ability to 
execute detailed vorticity measurements in such flows. The characteristics 
of many technologically interesting flows can be modelled in terms of the 
behavior of their mean vorticity , see e. g. , Foss and Kleis [1976]. The 
control volume form of the vorticity transport equations (Potter and 
Foss [ 1975 ], Chapter 5) shows the importance of the surface pressure 
gradient for such flows. The role of the time- dependent vorticity field 
in the production of acoustic noise has been developed in Eulerian form by 
Hardin [1973], This relationship expresses the far field density fluctuation 
p (^t) as a function of the vorticity and velocity of the flow field; 

Si 

specifically 



It is this last application which has been the principal motivation 
for the development of a capability to measure one component of the 
three-dimensional vorticity field; however, such a capability will also 
serve many other problem areas including those previously identified in 
this introductory section. 

Comprehensive and careful measurements of the mean velocity 

field would allow the mean vorticity field to be evaluated; this however 

would not meet the desired capahiKtyl An instantaneous evaluation is 

Z 2 * 

required if a quantity such as 9 ^ x V)/9t is to be determined. 

Consider a flow in which the gradient of the mean shear lies in 
the x-z plane. The transverse, or y- component, of vorticity is 

9u 9w 
^y ~ 9z 8x 

The measurement of this component is the object of our effort. It has 
been selected because it represents the principal component in the class 
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of flows under consideration {the mean values of w and w are zero). In 
particular, it is important in the construction of u x V for the flow field 
of particular interest, viz. , the normal impingement of an axisymmetric 
jet. It is shown in the following that this quantity can, in principle, be 
measured if one is willing to accept the approximation to a point measure- 
ment afforded by an array of probes and if one is willing to make further 
assumptions regarding the construction of the derivatives. The same 

probe configuration could be used to attempt a measurement of u , albeit 

2 

the approximations will be more severely tested by the presence of mean 
velocity gradients. The measurement of is much more involved in 
principle; however, a generally successful approximate technique is 
available; see Kistler [ 1952] or the report by Corrsin and Kistler [ 1955]. 

The essential strategy of the present measurement is demonstrated 
in the schematic representation of the vorticity probe shown in Figure 1. 

The quantity Wy. will be evaluated by the following procedure: (1) evaluate 

I V j and Y as the two unknowns which are related to and of the 
x-wire probe; (2) construct 8w/9t for the time (t+At) from the time series 
w(t) as 

— - w(t+2At) - w{t) . 

at - 2At ’ 

(3) construct 9w/ax (at t+At) by utilizing the "frozen flow assumption 

8w ^ 2. 8w 

8x u 8t 

(4) construct (8u/8z) from the readings of the parallel wires as 

(x, y, z) ^ cos y(x, y, z) {[V(x, y+Ay, z+Az) - V(x, y+Ay, z-Az]/2Az]} 

where the approximations resulting from the lateral displacement of the 
probe and the utilization of a single value of cos y are immediately 
obvious. 

The strategy of the measurement can be used to identify the 
specific considerations which represent the contributions of this report. 
First, a scheme to reliably compute | V[ and y from Ej^ and E^ must be 
developed. This scheme is the subject of Section 2. The operation of 
forming the difference [(8u/8z) - (8w/8x)] of two quantities which are 
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themselves differences suggests that questions of (1) accuracy and (2) 
uncertainty are of paramount importance in this study. If is to have 
integrity, or at least if the reliability of the final value 

is to be estimated, then considerable attention must be given to matters 
of accuracy and especially uncertainty. Section 3 presents a comprehensive 
evaluation of these matters. 

It should be noted that the preliminary considerations of the 
response equations and the uncertainty estimates, as given by the present 
author in the Second Semi-Annual Report (Foss [1975]), are obviated by 
the developments reported herein. Specifically, an improved scheme 
to convert (E^, E^) into (|v|,y ) has been developed and a considerably 
improved methodology for constructing the uncertainty considerations 
has been evolved. The earlier efforts were useful in the evolutionary 
process but they have been uniformly superseded by the present work. 

2. COMPUTATION ALGORITHM FOR u, w GIVEN 

The response of an individual hot-wire channel is dependent 
upon the magnitude and direction of the velocity vector with respect to 
the hot-wire itself. If the yaw angle is zero (that is, if the cross-product 
between the velocity and a vector parallel to the wire is perpendicular to 
the probe axis), then the response of the two wires constitutes two 
equations for two unknowns. It is the purpose of this section to develop 
the algorithms by which the hardware of the VORCOM or software in a 
general purpose computer can be utilized to extract the magnitude and 
direction of the velocity vector given the two voltage readings. 

2.1. Pitch Angle Response ^ 

Several pitch angle response equations have been offered in the 
general literature. A summary of the more prominent ones include: 

Hinze [1959] and Champagne et al. [1967] 

Veff = |v|[cos^a + k^ sin^ (8) 

Fujita and Kovasznay [1968] 

V^^£ = |v[[cos Q + e(cos a - cos 2a)] (9) 
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Brunn [1971] 

V^£f = Jv| cos“a (10) 

Friehe and Schwarz [196 8] 

^eff ^ lV| {l - b[l - a]}^ (1.1) 

where a is the angle between, the normal to the wire and the velocity 
vector. Note that each of these is based upon a modification of the 
"cosine law" which would be valid for a wire of infinite aspect ratio . 

For the purposes of the present study, it was considered 
satisfactory to identify an effective pitch angle response formulation 
which fit the empirical data and was readily integrated into our computation 
scheme; the Friehe and Schwarz relationship was selected in response to 
these criteria. 

E. 2. Solution for the Pitch Angle 

The computational strategy to recover a given Ej and E^ will be 
established in this section. The pertinent terms of the equations will 
first be defined. 

A two-step data acquisition process will be utilized. The probe 
will first be aligned with the time mean flow; a time series of u, w, and 
Wy 'will then be obtained. For an axisymmetric flow, the probe will be 
positioned in the x-r plane. The instantaneous velocity for the time 
series data will-be described in terms of its pitch angle (y) with respect 
to the probe. (The final data will, of course, be referenced to the x, r, 

6 laboratory coordinates.) Consequently, as shown in Figure 2, the 
pitch angles with respect to the hot-wires of the "x" array caii be 
expressed as 

Ui = -Y and = p^ + Y (12) 

where the p values are defined on the figure. The effective velocities are 
related to the hot-wire voltage values by the relationship 

+ K. (o) }=1,Z (13) 
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where and K are constants for a given overheat ratio and ambient 
j 

temperature. A large number of calibration data sets generated in our 
laboratory have supported the original proposal of Collis and Williams 
[ 1959 ] that m = 0. 45; however, the most appropriate value of the constant 
m is determined for each data set by the calibration process described 
elsewhere in this report. 

The essence of the computational strategy is to compute y from 
the voltage pair (E^jE^) and equations (11), (12), and (13). An explicit 
form for the relationship between Ej^, E^ and y can be developed as 

(E^ - |{l-b2[l-cos^/^(P2-Y)]} J 

2f 

by using (11) in (13) and cancelling | Vj . 

The left side of (14) will be referred to as G and the right 

side as <G> indicating that the right side is computed from the results of 
the prior averaging processes which yield K(y)* It should be noted that 
b, = b^ = 0* 92 is assumed on the basis of the Friehe and Schwarz data 
and that = 0. 45 will be used unless a particular probe calibration 

appears to require different values, pj and p^ will be individually evaluated 
from the symmetry_point of the hot-wire angular response curves. 

The magnitude of the velocity can be determined once the pitch 
angle is evaluated. Specifically 

|vl = [Ej - E^ * F(y) (15) 

1 

where F(y) is defined as 

F(y) = K(v)"^/“^ {l-b [ - cos^/^ (p-V)]}~^ (16) 

If the G value was everywhere equal to the <G> value, then 

meas 

the above considerations would be sufficient for the calculation of y and 
j v| . However, the calibration of Ej(| V|, y) and E 2 (l V], y) require the 
definition of the following quantity 
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G 


meas 


- <G> = A(v) + S(y) (1 V| - <V>) 


(17) 


where the reference velocity of <V> is arbitrarily chosen as an approximate 
center point of the calibration data set. An iterative scheme for the 
assessment of the y value is therefore suggested. The details of this 
computation and the calculation steps to be executed by the VORCOM are 
shown in Figure 3b, The schematic diagram showing the steps in the 
computation process is considered to be self-explanatory; the ability of 
a single iterative loop to evaluate the correct y is discussed in the 
following section. 

2.3. Characteristic Results 

Several complete data sets have been obtained: Ej^(|v],y), £^(1 Vj,y), 
for 10 - |vj 2 120 where 

!^fi+i = l^li + 

and for - 40 ^ y ^ 40 degrees where y^^^^ = y^+ 5 degrees are representative 
of the data base required for a complete set. These numbers have been 
used to establish the validity of the basic form of the response equation 
(including the determination of E^) and to evaluate the functional form 
for K(y) from which the following are obtained: 

<G(y)>, F(y), S(y), A(y) . 

Graphical representations of these relationships are presented in Figures 4 

through 6 and Table 1. The A(y) and S{y) values of Table 1 suggest that 

these quantities are not a smooth function of y. This lack of 

smoothness is attributed to their small numerical values and hence to 

their susceptibility to inaccuracies; however, this is not deemed to 

represent a significant degradation in their utility for the calculation 

scheme. Table 1 also includes the standard deviation (S.D. ) of the 

quantity GDIFF where GDIFF = G - <G> and 

meas 

N 

^ actual ' A(v) + S(v) (V. - 

i=l 


(Note that <V> = 55 fps is used for this data set. ) 
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This standard deviation is, in some cases, a significant fraction of GDIFF 
itself; however, in all cases the quantity which represents the effect on 
the calculation of y is small. Specifically, iet represent the 

error resulting from the use of a linear approximation of GDIFF, then 


error 


- ^ 
dG 


S.D. 


The values of y are also presented in Table 1; these small values 

error 

justify the adoption, of the indicated calculation procedures. 

The G values are identified in Figure 5 for the set of 

meas 

velocities 10 ^ lv| ^ 120 fps. From this presentation, it is apparent 

that G is quite close to <G> for the small absolute values of y and 

meas ^ 

that the deviations between G and <G> are small but not a smooth 

meas 

function of y. The linear form used in (17) was selected for computational 
convenience and because of the monotonic dependence of the difference 
in the G values with respect to V. The pertinent feature is the ability to 
recover the correct y value from the and E^ values. An evaluation 
of this "recovery ability" is included in the discussion of "accuracy, " 
see Section 3.1. 

The ability to recover ] V| given y is a more precise process. The 
validity of theK{y) evaluation from the initial calibration data sets is felt 
to be essentially limited by the accuracy of the V calculation based upon 
the measured total pressure. Since Kis a smooth function of y, and since 
the resulting F(y) is also smooth, there seems little additional uncertainty 
in these calculations once y is known. This assertion is also evaluated 
in Section 3. 1. 


3. ACCURACY AND UNCERTAINTY 

The central question to be addressed by this section is: 

Given E^, E^,, E^, E^ as voltages from the hot-wire 
channels and given the calibration data and computational 
procedures described in the previous section, what is 
the relationship between the calculated u, w, and to^ 
values and the values of the same quantities which exist 
in the jet flow. 

The response to this question is best constructed in terms of the 
accuracy and the uncertainty of the measures. It is useful to first provide 
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an explicit definition of the two words. ^ 

accuracy: degree of conformity of a measure to a standard 

or true value 

uncertainty: not having certain knowledge 
We will here speak of the accuracy of the computation of u, w, lo given 
^ 2 ^ ^2^ ^3 '^known^' conditions of the flow field as the accuracy 

of the computing procedures. There are three aspects to the designation 
of the accuracy. First, the calibration data, which are used to compute 
K(*y) s-^id p, are known at discrete values of the velocity and the pitch 
angle y; hence, it is necessary to construct interpolation schemes for the 
construction of intermediate values. Second, the control variable V is 
not precisely known during the calibration process; that is, the total 
pressure required for the calculation of V is limited by the available 
accuracy of the capacitive pressure transducer. Third, the accuracy of 
the hot-wire voltage reading is limited by the A/D converter resolution; 
specifically, the voltage span of 0 < E 4 volts will be processed by the 
12 bit A/D converted of the Texas Instruments minicomputer (in the 
calibration process) and the approximately 1-1.2 volt active portion of 
the signal (E^ < E < will be processed by a high speed (4 psec) 

10 bit A/D converter for the data acquisition. In both situations, the 
expected resolution is of the order 1-1.5 mv. Significantly, the noise 
of the hot-wire anemometer is < 1.5 mv; hence, noise will only effect 
the least significant bit of the reading* It should be noted that this noise 
level includes the beneficial effect of the 20 \isec averaging time; this is 
essentially a low pass filtering operation and, as such, it removes the 
high frequency noise of the anemometer circuitry. Spurious effects such 
as d. c. drift of the amplifiers, changes in the ambient temperature level 
and dirt accumulation on the probes will be limited by careful monitoring 
of the experiment. The specific aspects of the accuracy of the computing 
process are evluated in subsection 3. 1. 

When the four-wire probe is located in the jet flow, the factors 
which result in the observed voltages cannot be uniquely identified in 


Webster s New Collegiate Dictionary 
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terms of the values V and y. That is, a transverse velocity component 
V may be present in addition to the u and w components in the plane of the 
x-wire. Secondly, the possibility of a gradient in the pitch angle y and/ or 
the velocity magnitude V must be recognized. These effects cannot be 
ascertained using the signals of the x-wire array. However, the effects of 
these additional factors can be analytically modeled and their magnitudes 
can be represented by appropriate statistical measures. Because the 
instantaneous values, which will comprise the time series for u, w, and 
Uy, are influenced by effects which are unknown, i. e. by factors which 
are represented by our ”... not having certain knowledge . . . , ” we refer 
to these effects as representing an uncertainty in the values for u, w, and 
Wy. These considerations are discussed in full in subsection 3. Z. 

3.1. Accuracy 

A sequence of operations is required to assess the vd.ocity in the 
flow field. The essential feature which defines the accuracy of this process 
is the identification of the least accurate step in the sequence. The three 
major steps in the sequence are identified in the opening discussion of 
Section 3; they are (1) reliability of the interpolation formulae, (2) measure- 
ment of velocity in the calibration process, (3) A/D converter resolution 
and the noise level of the hot-wire anemometer. 

The measurement of the velocity in the calibration process is considered 
to be the limiting factor in the accuracy of the measurements to be made 
with the VORCOM. Specifically, the linearity specification for the 
capacitive pressure transducer and the nonlinear V ~ N/p relationship 
results in a nearly constant 0, 5 fps resolution in the velocity. This 
a priori assessment has been indirectly verified by numerous calibrations 
of various wires and the subsequent calculations of the standard deviation 
of • - V ] . A typical set of calibration data is reproduced in 

Table 2 and is presented graphically in Figure 7. 

The accuracy of the interpolation relationships described in Section 
2 is difficult, to assess given the recognized accuracy limits of the original 
velocity measurements. However, it has been found that the functional 
form 

= E^ + K (18) 

o eff 
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fits the calibration data to within a nominal value (standard deviation) of 
~0. 6 fps if a constant value of m is used for the fit; see Table 2. The 
calibration data invite the following rather interesting speculative 
considerations regarding the recoverable accuracy of velocity measure- 
ments, Specifically, if the above relationship (18) is accepted as valid, 
and if the hot-wire voltage measurement is accurate (to within the noise 
it at ion of the instrument) then the hot-wire voltage can provide velocity 
measurements with an accuracy exceeding that of the original calibration. 

That is, the standard deviation quoted above reflects the uncertainty in 
the velocity as monitored by the pressure transducer, a more accurate 
measurement of would conform to the deduced from the 

^o values of the calibration process. 

The precision of the 10 bit A/D converter and table-look-up operation 
of the VORCOM is decidedly not the limiting factor in the sequence of 
operations. Specifically, the conversion/look-up process can provide 
— 0, 1 percent resolution over the full range of values; the accuracy of the 
measurement is given by the estimated :t 0. 5 fps ambiguity in the calibration 
velocities which corresponds to ± 5 percent accuracy at 10 fps and 
± ~ 0, 5 percent accuracy at 120 fps. Hence, the accuracy of the measure- 
ment is limited by the accuracy of the original velocity measurement and 
not by the resolution of the signal processing equipment.’ 

Since these estimates of the accuracy are based upon a number of 

assumptions, it is useful to also characterize the accuracy for a given data 

set in terms of the ability to recover the known | Vj and (y). This comparison 

is presented in Table 3. The calculated values are derived from the 

and E^ values of the calibration data set; they are respectively compared 

with the y, which was set in the calibration process, and the V, which was 

calculated from the pressure transducer. As noted above. Section 3.1, 

deviations up to ± 0# 5 fps can be expected as a result of the pressure 

transducer characteristics. The comparison between |v| , and |v| 

^ ' 'calc ' 'meas 

and that between ^nd presented as both a percentage and 

as an absolute value for each velocity, pitch angle pair. 


Brunn [1971] makes use of calibration data over the range 0. 4 ^ V ^ 1 50 mps 
to evaluate m(V). In the range of current interest 3 ^ V ^ 40 mps, a 
a constant value of m is also supported by his data. 
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The accuracy level implied by the ± 0. 5 fps estimate is quite 
adequate for the determination of the velocity magnitude. As will be 
evident in the following section, the inescapable uncertainties associated 
with the measurement of u, w, and u in the hostile environment of a 

' ' y 

turbiilent flow will be the limiting factor in the interpretation of the 
original data. In this regard, it is especially important to note that the 
differencing operations of the voriticity evaluation are sensitive to whether 
the accuracy limitation results from erratic or smooth perturbing effects. 
The nonlinearity of the pressure transducer and the inaccuracies of the 
interpolation formulae are all "smooth" in character. Hence, the 
contribution of the inaccuracy will tend to cancel out with the differencing 
operation. Electronic noise would, of course, be erratic and would make 
a statistically significant contribution were it not for its small magnitude 
(= 0. 1 percent of the voltage measurement). 

3.2. Uncertainty 

The factors influencing the uncertainty will be traced by following 
the signal processing sequence to yield the u, w, and values {referenced 
to the probe coordinates). The following convention is introduced for this 
analysis: 

i. The measured value of the quantity of interest, [ ], will 
be expressed as its 'true" value, [ ].p, plus a difference 
term, 6[ ]. 

[ ] = Ux + 6[] (19) 

ii. The relative value of the difference term is expressed as 
e|- j and is defined by the expression 

€[ ] = [ ]/[ ]-p - 1 • (20) 

In the context of this discussion, the "true" value, [ ]^, will not 
include the considerations of accuracy as presented in Section 3. 1. That 
is, we seek only to describe the effects of the actual flow field which 
result in our "not having certain knowledge" as regards the relationship 
of the measured to the true value of the quantity [ ] . 
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The meaning of .the uncertainty can be further defined by noting 
that (i) the uncertainty of the calibration data is zero and (ii) the causal 
factors of the uncertainty are known but their magnitude and hence their 
influence on each individual reading are both unknown and, for the 
experimental capability at hand, unknowable ♦ 

The strategy for the evaluation of the uncertainties follows from 
the second condition; namely, the (uiiknowable) effects in the flow field 
which influence the magnitude of the measured quantity [ ] will be 
analytically incorporated into the calculation formulae for [ ]. The 
magnitude of the disturbing effect will then be characterized by its 
standard deviation to allow the value of Cj- j to be evaluated. The final 

results will be presented as ey and where the j values 

are functions of y ^^d are parametrically dependent upon the quantities 
Ay, A V/V and v/V which are defined below. The values are dependent 
upon the results of special experiments to be executed; however, estimates 
of that fraction of the values that can be extracted from these purely 
analytical considerations are also presented in the following. 

The perturbing effects, which create the uncertainties, are a 
result of the three-dimensional, spatially nonuniform conditions of the 
actual (turbulent) flow fields. Specifically, the magnitude of the velocity, 

[ V|, and the pitch angle, y, may be different at the two wires of the 
x-probe. In addition, it can be expected that a lateral velocity, v, will 
be present in the flow and hence | V| = [u^ + as compared 

with the calibration condition in which | V| = [u^ + 

These three effects will be characterized by the quantities 

3,nd (v/V) where the A quantities are the difference values 
between the two wires. 

3. 2, 1. Uncertainty in y 

The pitch angle y is dete r mined from the magnitude of the G 
function. The uncertainty in y can be determined by evaluating the effect 
of the three perturbation quantities. Ay, (AV/V) and (v/V), on the magnitude 
of G and hence on the magnitude of y. The latter relationship can be 
expressed as 6y (where 6y = y^ - y) 
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and 


6y = Ig SG 


(21) 


6G = G - (22) 

The quantity G represents the value obtained by the computation procedure 
described in Section 2 ; that is, Yj F(Yj) '*’ ! ^1 ^ l-*^meas “ 

-»■ <G(-Y 2 )^ “*■ Y 2 ■^G(‘Y 2 )^ is the G value of equation (22).* The true 

value of G (i. e. , G,p) is (defined to be) that value which would have been 
produced by the same steps if Ay, (AV/V) and (v/V) were identically 
zero. Since an explicit relationship for G is available, the quantity 6G 
can be computed as 

6G[y, 1v|;Ay,(AV/V), (v/V)] = G[y, | Vj; Ay, A V/V, v/V] 

-G[y, |V|] (23) 


where 

G[Y.|vi;Ay, (AV/V), (v/V)] = 

Ki(Yi)^/*”l[{l-b[l - co s^^^(P^-Y^)]}'^ + (v/V)^]^^^ 

K2{Vi+Av)‘^’"2[{1-K1-cos'^^(P2-Vi-^Y)]}'^+{vi+Av/(Vj+&V)}^]‘/‘‘(Vj+AV) 

(24) 

An immediate simplification appears to be in order, viz. , Av = 0 will 
be assumed. 

The desired quantity, 6y, follows from equations (21), (22), (23), 
and (24) where G[y, | V|] is obtained by setting Ay = AV = v = 0 in (24) 
and where the explicit dependence upon V can be eliminated as seen by 
(24). 

3. 2. 2. Uncertainty in V, u, and w 

The calculation of V requires a known value for y and the voltage 
The computing equation is based upon the assumption that the velocity 
is in the x-z plane; that is, that the conditions of the calibration process 

Note that the subscripts 1 and 2 refer to the initial and first iteration values 
for y; see Figure 3b. However, the subscripts in (24) refer to wire 1 and 
wire 2 of the x-array. 
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are duplicated in the flow field. Hence v exerts an implicit (6 y = 6y(v/V, . . . )) 
and an explicit perturbing effect on the [ V | uncertainty . The development 
of Cy can be constructed by approximating with the expression 

( 25 ) 

presumes tha.t the response equation (18) is unequivocally valid and 
that yrj, is the pitch angle at wire 1. The unequivocal valididty of (18) 
implies that the pitch angle relationship, including the value of K, is not 
effected by the presence of a non- zero yaw angle. Hence V„ of (25) is the 
velocity magnitude in the plane of the x-wire, viz. | | = (u^ + w^)^/^. 

The value of | V| which is determined by the computing equations 
is based upon the measured E^, the spurious y (i. e. , + 6y) and the 

neglect of the lateral velocity v. This value can also be expressed using 
equation (18), viz.. 


E^ = +K(y^ + 6y){l-b[l - cos^/^P^-y-6Y)]}^"^|v|"^ (26) , 

2 

By subtracting E^ from both sides and making use of the definition of €|- j, 
see (19), the quantity Cy can be expressed as Cy = V/V„ - 1, hence 


1 + e 


V 


~K(Vt) [{l-btl-cos^/^(p^-V^)]}^^. (v/V^)^]*/^ 

K(Yj+6y)J {l-b[l.cos*/^(Pj-YT-6Y)]}^ 


(27) 


The uncertainty estimates for u and w follow from the Cy and 6y 

values which have been established above. The expressions for e and e 

u w 

are based upon the definitions for u and w, viz. , 

u = V cos Y and w = V sin y (28) 

Substituting u^ + 6u for u, V.^, + 6V for V and y.j. + 6y for y, yields 
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+ 6u = (V^ + 6V) cos(y^ + 6y) 

= (V^ + 6V) (cos -y^ cos 6 y - sin y-j. sin 6 y) 

= cos Y.J, -V^y sin y^ + 6V cos y^ - 6y sin y-j 

= cos Y.J. [ 1 - 6 y tan y^ + €y - €y 6y tan y,j. ] 

or 

= [cy - 6YtanY^ - eySYtanY^] (29) 

Similarly, 

W.J, + 6w = (V + 6V) sin (y,j. + Sy) 

= (V + 6V) [ sin Y-p cos Sy + cos y^ sin 6y] 

= V sin Y^ + V 6y cos y^ + 6V sin y.^, + 5V 6y cos y^ 

= V sin Y^ [ 1 + 6y cot y^ + Cy + Cy 6y cot y^ ] 

or 

[ Sycot Y^ + €y + €y Sycot Y^] (30) 

Representative conditions for (Av/V), Ay, and (Av/ V) have been used to 
compute numerical values of t he uncertainties 6y, e^, and e^. The 
results of these calculations are summarized in Figure 8. 

3.2.3. Uncertainty in Uy 

The vorticity component transverse to the time mean streamline 
(and hence to the probe axis), , is defined as 

3u Ow 
^y ~ 9z 9x 

The y. component vorticity is constructed as the difference of two quantities 
which are themselves differences; numerous factors must be considered 
in its evaluation. 

The uncertainty, 5i0y, is defined by the expression 

Wy = Wy]p + 6w^ (31) 
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Since the uncertainty can also be expressed as a Taylor expression about 
the true value w L , the form 

y 




“y^T 


+ 


8c0y 

9(8u/9z) 


a 9 w 

+ . y , 

^8z'' ^ 8(8w/8x)- 


(— ) + 
'8x ^ 


( 32 ) 


is appropriate. Hence, 


8c0y = 6(8u/8z) - 6(8w/8x) 


( 33 ) 


if the linearization of the Taylor series is valid, that is, if 6 ( 9 u/ 8 z) and 
6 [ 9 w/ 9 x] are sufficiently small. Such a condition will be assumed. 

It will be convenient to further subdivide this subsection in order to 
allow separate considerations of the uncertainty in 8w/8x and 8u/8z. 


3 . 2 . 3 . 1 . Uncertainty in 8w/8x 

The quantity 8w/8x can be formed from the instantaneous z-component 
Navier-Stokes equations. Specifically, 






8w , 8w , 8w 
-sr~ + V -5 — + w -Q — 
8t 8y 9z 



8z 


8^w 1 

8x. 8x. I 
J J-* 


( 34 ) • 


If the "frozen flow" assumption, which is implicit in the Taylor hypothesis, 
is made then the following formulation is a rational construction for 
8w/ 9 x 


8w _ 8w 

8x “ ~ u 8t 


( 35 ) 


An expression for the uncertainty will be extracted from ( 34 ) and 

( 35 ); defined by the expression • 

If = 

Substituting ( 34 ) and ( 35 ) into ( 36 ) yields. 


u 8t 


1 

u 


9 w , 8w 

ar + '^af 


+ w 


8w 

9 z 



9 z 


8 ^w 

9 x. 9 x, 
J J 


^ 9 / 9 x^ 
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or 


1 = 


1 + 


8w , 9w , 1 9p 
9y 9z p 9z 


9^w 


9w 

9t 


9x. 9x. 

__! L 


^a/9x^ 


(37) 


Considering that Cgy'gjj 1 following approximation for Cgy'g, 

■fae derived 


can 




9y 


(38) 


1 J 


The expression for 9w/9x shows that this quantity is also influenced by 
the uncertainties associated with u and w. Consequently, the uncertainty 
^9w/9x defined by the expression 


9w 




9x u^(l+e^) 9t 


+ ^9/9x^- 


The quantity 9wrp(l+e )/ 9t can be operated upon to yield 
i. . W’ 


(39) 


at " at 


{l +6 H" 




w 8wrp/ at ' at 


Be 

(^)> 


(40) 






9w„ 


9t 


9t 


■ 9w„ 
(9t ) 

9t 



w 

T 


96w 

6w 

9w— 

9t 

2 9t 

<8t ) 


Wt 

96w 




at 

6w 



9w^ 

Wt 




= (l^wl-^w) 


(41) 


where the last step involves the central assumption that the ratio of the 
time derivatives is adequately represented by the ratio of the arguments 
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of the differentiated- quantities, namely 6 w/w_ = [e |. Therefore, 


w' 


8w_(l+e ) 8w 


at 


■w'" at 


( 42 ) 


and 


nj^ \ 

dx ' ^dw/d^^ 


X 


1 

J i ll + 

at ^ 




_L ^ 

at 


} 


{l +[e + € + € 

^ w u 


9/dx- 

2 jl/2y 

(43) 


a/8: 


It should be noted that will have to be evaluated from supplementary 

experiments as described in the following. 


The evaluation of the quantity e 


3/d: 


required the evaluation of^ 


spatial derivatives of velocity and pressure, see equation (38)* For 

convenience, the spatial derivatives will be grouped into three quantities 

to be separately discussed; the groups are (v 3w/ay + w3w/dz), p"^ ap/3z, 

and i/d w/dx. dx., 

J J ■ 

The magnitude of the kinematic viscosity times the sum of the 
2 

second derivatives (i/ 3 w/dx. dx.) is expected to be small; specifically, 

it iS' expected that the second derivatives themselves are small. The y 

2 2 . 

and z derivatives are not accessible,* the quantity 8 w/dx can be 

approximated using the same Taylor hypothesis which is presently under' 

investigation* If it is assumed that the' small scales of the motion which 

2 2 

are responsible for 9 w/9x are essentially isotropic, then 


a^w 


dx. dx. 
J J 


(2^^2' + 1) = 3. 

. ax^ 


83 


a^w 

ax" 


(44) 


and, from the frozen flow approximation. 


A total of six channels of anemometry and three x-wires would be required, 

i. e. , a w/3y = [w(y+Ay) - w(y-Ay)]/2Ay would require two x-wires and 
2 2 

four channels, d w/dz could be obtained from ah additional x-wire at 
(y, z-fAz). 
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(45) 


9^w 


8x 


3.8[^ ^ 

u at 


1 dw 3u 1 

3 at at 

u 


The mean square pressure gradient d-p/dz in isotropic turbulence 

has been analytically related to the longitudinal correlation fimction, see 

Hinze [1975], p. 308. Since the pressure gradient is sensitive to the 

small scales of the motion and since, as in the above assumption for 
2 2 

a w/8x , local isotropy is (perhaps) reasonable. We will make use of 

2 

this relationship to estimate (8p/az) 


(92) 

^ 3z^ 


= (^)' 
'8x^ 


= 4 


Z I 4 
P ^ 


I 


CO 


f(r') dr' 


(46) 


where u is the r.m. s. value of the longitudinal velocity and f(r) is the 
longitudinal correlation function. By again invoking the frozen flow 
hypothesis for the construction of f(r), a time series for u(t) is sufficient 
to evaluate the mean square pressure gradient. 

It is expected that the quantity (v 8w/ 8y + w 8w/ 8z) may be 
quantitatively significant with respect to dw/Bt; hence its magnitude must 
be assessed as accurately as possible. A complete measurement of this 
quantity by finite difference techniques is not feasible since three x-wire 
probes would be required for v8w/ay alone and this measurement should 
be made simultaneously with w8w/8z which would require two additional 
x-wire probes. . . 10 channels of anemometry and a "forest" of hot-wires. 
Consequently, the experimental evaluation of the desired quantity must 
necessarily involve rather substantial approximations. Several alternative 
schemes were considered; the following is deemed to represent the optimal 
evaluation of the quantity. 

Consider that an x-probe is oriented to measure the z component 
velocity and that it is placed at (x, y + Ay, z + Az); this location will be 
referenced as +A, A second x-probe is located at (x, y - Ay, z - Az); 
this location is term -A. The z component velocity at +A is related to 
that at -A by the expression 


w(+A) 


w(-A) + EAy + 


9w 

8z 


2Az + 


a w 

a 2 

8y 


(2Ay)‘ 

21 


a w 

8z" 


(2Az)^ 

21 


, 3w 
8y aw' 


AyAz + • • • 


(47) 
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The quantity 8w/9y + 9w/3z can be extracted as 


^ w(+A) ~ w(-A) 

9y 9z 2Az 


( 48 ) 


where Az = Ay is ensured by the placement of the probes and where the 
higher order terms R will be as sumed to be negligible. The assumption 
is rational if the probe displacement is sufficiently small. The desired 
quantity differs from that given in (48); viz.. 


V 


9w 

9y 


+ w 


9w 

dz 


versus 


9w 9w 
9y 9z 


The above epxerimental arrangemert allows w(x, y, z) to be approximated 
as 


■w(x, y, z) 


w(-fj^) -)- w(— 


(49) 


and, if the further as sumption is made that w ~ v, then 


9w 

^8y 


+ w 


9w 

8z 


f . ,0w , 9w. 
W(x, y, X) 


^ w(-hA) + w(-A) .w(+A) ^ w(-A) , 
2 2Az 

- - w^(-A) 


(50) 


The assumption that v - w is conservative since it is known that vw = 0 . 
Hence, the uncertainty €g can be written as (see (43)) 




X 


= {fv-^ + w— + [y — ^ + r- 

^L^ay ^^az-* ^ LI' 9x. 9x. •* az-* 

J J 


_ / . p w (-}A) - w (-A) 1 ,3.8 ,w(t+At) + w(t-At) - 2 w(t) ^ 

'ZAz' *• w(t+At) - w(t-At)-* ■ ^ L w(t+At) - w(t-At) ' 

. ,3. 8 V 2 „(t+At) - , 4i^4"[f(r')/r’] dr' 

+ [ —) L 2St J ^ 


{[w(t+At) - w(t-At)]/2At} 

(38a) 
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wTiere the assumed statistical independence of the three groups of terms 

2 

requires that they contribute to (^ 0 ^ 0 ^) as the sum of their squares. A 
review of the construction of ^ 0 ^ 0 ^ suggests that the several major 
approximations should lead to a conservative estimate for the uncertainty. 


The full value of the uncertainty in e 


8w/ 9 


X 


can now be constructed. 


Specifically, the relationship for from (38a) can be utilized in the 

expression for from (43). This lengthy expression will not be 

explicity written here. 

An examination of equation (43) reveals that a portion of it may be 

2 2 1/2 

calculated a priori; viz., [e + e ] ' . This quantity will be termed 


W 


it is included in Figure 8 for reference purposes. 

3. 2. 3 . 2. Uncertainty in 9u/9z 

The quantity 8u/9z will be constructed from two single wire probes 
which are parallel to the y axis, displaced a distance 2 Az apart, and 
located at a distance of Ay from the x-wire. The x-wire occupies the 
location of interest, viz. , (x, y, z). The essential strategy of this 
measurement is to determine the desired quantity dy/dz at x, y, z from the 
value 


V(x, ytAy, - V(x, y+Ay, z-Az) 


(51) 


where the single wire is assumed to respond to the magnitude of the x-z 
plane velocity component. (Note: This assumption of the cosine relation- 
ship for the yaw (i. e. , v) effects of the straight wires is justified on the 
basis that it is less restrictive than other assumptions which will 
necessarily be made. ) 

The measured quantity is an admixture of information from the y 
and the y + Ay planes; consequently, one element of the uncertainty is the 
magnitude of the ratio 


V (x, y+Ay, z+Az) - V (x, y+Ay, z-Az) 
V(x, y, z+Az) - V(x, y, z-Az) 


(52) 


If the central difference value were an accurate method of constructing 
the derivative and if there were no uncertainty in assessing or utilizing 
cos Y, then 


9u 

9z 


9z^t ^ ^ 
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( 53 ) 



would relax to 


9^ - n 4. 

^ - d^L ^ 


Ay 


] 


(53a) 


Note that the numerator of (52) is 'the measured value whereas the 
denominator is equivalent to the true value for the stated conditions. Hence, 
the quantity (1 + is equal to the ratio given, in (52), A rather simple 
procedure exists for the evaluation of this ratio; a four wire array with 
single wires parallel to the y axis will be made and the r.m. s. value of 
the ratio will be used to compute Namely, 



V (x, y+Ay, z+Az) - V (x, y+Ay, z-Az) 

2 1 

1^0 Jo 

V(x, y, z+Az) - V(x, y, z-Az) 

d 


(54) 


The motivation for using -y(x, y, z) versus ^(x, y, z+Az) and*y(x, y, 
z-Az) is simply that the first value is available during the measurement 
of Wy . An expression for the uncertainty associated with this aspect 
of the measurement can be developed as follows. Since u = V cos y 


9u -1 

J 


9z ■'x, y, z 9z 


= CO® y] 


x,y,z 


= cos y -^1 + V » 

’ 9z ■'x, y, z 9z 


9 cos ^ j 


X, y, z . 


— CO s Y "q — ] ,A 

' 9.Z -'x', y+Ay, z ’• Ay^ 


1^] (55) 

oz y, z 


where, the first term on the r.h. s. contains the measured quantity and 
the formulation 


9u- i _ 9u IP 1 ■ 

9z-'measured ” 9z-’T^ ^ ^A cos • J (53b) 

and a manipulation of (55)* allows the following definition of e. 

Acos 


Note that d\i/dz)j, is the 9u/9z on the lefthand side of (55). 
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cos y) 


X, Y, a 


9V . ^ ^1 

9z ^x, y+Ay, z 8z -^x, y, z 


1 


.. 9 cos Y ^ 

9z 'x, y, z 

9u > 

9z -^x, y, z 


1 + e 


AyJ 


or 


V 9 cos y/9z) 


cos 


9u/9z) 


X, V, z 


X, y, z 


(56) 


An estimate of e can be extracted from the data base used to 
cos 

evaluate the convective acceleration terms of the Taylor hypothesis. 
Namely, from two x-wire probes measuring u, w at (x, y-Ay, z-Az) and 
(x, y+Ay, z+Az) respectively, we can make the following approximations 
(recall that +A implies x, y+Ay, z+Az and -A implies x, y-Ay, z-Az): 


u(+A) - u(-A) _ , 9u 9u . 

2 Az ^ 9y 9.z^ 


(57) 


and 


cos y(+A) - cos y(~A) _ 9 


2 Az 


9y 


cos Y + cos Y 
‘ 9z ' 


(58) 


By squaring and time averaging each side of (57) and (58), we obtain 

< = < [ (is)' + <ls> " 


9y' 


(59) 


and similarly 


^|- y(±A)^_.V;(-A) cos V(+^^- V ^ ~ 2 (60) 


9z 


where the rather major assumptions of local isotropy allows the 
approximations 


.9u . 2 

'9y^ 


(— )^ 
'9z^ 


and 


'l7> 'H> = 


(61) 
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a.nd similarly for the (V cos y) product. Note that <[ ]> is used inter- 
changeably with [ ] for the time average of [ ] . Consequently, the output 
of the two x-wire arrays can be used to evaluate Specifically, 


e 


cos 



V(+A) + V(-A) 
2 


n 2 


cos y(4-A) - cos y('-A) 
u(+A) - u(-A) 



( 62 ) 


It is. somewhat discouraging to note that [cos y (+A) - cos y(-A)] will have 
to be quite small in order to overcome the large multiplicative factor in 
the balance of e 

cos 

The cosine term itself also introduces some uncertainty; viz. , 


cos y = cos (y,^ + 6y) 

and 

cos y = cos y^ cos 6y - sin y^ sin 6y 
= cos y^ - 6y sin y^ 

= cos y^ [ 1 - 6y tan y^ ] 


or 


e = 6y tan y 

cos y • ' 


(63) 


and ^ can be calculated using the previously identified uncertainty 

parameters: Ay, AV/V and v/V. 

The separate measurements and/ or calculations of the quantities 

^Ay’ ^cos’ ^cos y that they be combined as if they were statistically 

independent. Hence 


3u. X 9V, 

8^) = cosy)^) 


9z 'x, y+Ay, z 9 


’^T ^cos ^cosy^f 1 + e ^ 


Ay 


9u 

az->T 


Jrr, [1 + e 


8u/ 8: 


which serves to define the quantity e 


9u/ 9 s 
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3. 2. 3. 3. Statement of Uncertainty for 

The separate evaluations of 3.nd ^g^ 9 y can be collected 

to define e . Specifically, 

8u 9w 
‘^y ■ dz ~ 8x 


0z^T ^ ■ 8x ® ^8x^ 


8u . 8w . 

8z^T " 8x^T 


m|S) 




VtI ' + 


—\ ^ 
8z^T 8z 

‘^y^T 


8u ^ 8W '^ 8u > y8w . 

8z'T ■ 8x^T 8z^T “ '8x^T 

8w . c 8w '\ 

8x'T 8x I 




^y>T 


It will be consistent with the approximate nature of the e quantities 


to construct e^j using the (constant) values of ^8w/8x 

evaluated from^&e above described experimental data and calculations 
and multiply these by the measured 9u/8z and 8w/8x values. 
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cast epoxy 
probe body 



c 



Notes; a 
b 
c 


= .0. 2 mm 
= 0.4 mm 
= 2. 0 mm 


a - 45 degrees 

active portion of wire 5 [x dia. , 1 mm length 
total length of wire ~ 3 mm 


Probe- design based upon recommendations of Strohl and Comte-Bellot 
[1973]. 


Figure 2. Probe characteristics and definition of angles. 
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Figure 3b# The VORCOM* 
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Figure 7a. E versus V, hot-wire calibration. 

Note: Straight lines drawn through the discrete data points of Table 2 
constitute the "curve" which is shown. 
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Figure 7b. E“ versus m = 0. 45 

Note; Straight line segments have been drawn through the processed 
discrete data points of Table 2. 















Table 1. ' "Constants" for the y Computation Given and Values 


GAMMA 

A 

S 

SD 

(DGAM/0G)*SD 

-40 

1.26E-2 

-8.47E-5 

1.07E-2 

5.18E-1 

-35 

-2.51E-2 

-1.63E-5 

6.81E-3 

7.34E-1 

-30 

-2.98E-2 

9.27E-5 

2.52E-3 

2.00E-1 

-25 

-1.84E-2 

1.07E-4 

2.80E-3 

2.16E-1 

-20 

-4.53E-3 

9.42E-5 

9.32E-3 

6.24E-1 

-15 

7.03E-3 

l.llE-4 

4.44E-3 

2.44E-1 

-10 

1.48E-2 

6.62E-5 

2.70E-3 

1.20E-1 

- 5 

1.57E-2 

-1.16E-4 

1.59E-2 

5.68E-1 

0 

1.03E-2 

-2.11E-4 

8.49E-3 

2 *46E— 1 

5 

-2.16E-2 

-8.39E-5 

U02E-2 

2.42E-1 

10 

-3.99E-2 

-1.99E-4 

1.12Er2 

2.18E-1 

15 

-5.67E-2 

-3.03E-4 

8.07E-3 

1.29E-1 

20 

-3.10E-2 

-3.49E-4 

9.69E-3 

1.26E-1 

25 

2.26E-2 

8.56E-5 

2.27E-2 

2.29E-1 

30 

1.23E-1 

1.61E-3 

7.09E-2 

5.08E-1 

35 

1,48E-1 

5.54E-3 

1.95 E-l 

8.05E-1 

40 

-8.54E-1 

1.55E-2 

3.60E.-1 

-5,19E-1 


Notes:' (1) A and S are defined by the ejqiression: 

<G(y)> - = A(y) + S(Y)[V - < V>] ; 

see equation (14) for the definition of G. 

(2) SD = standard deviation between actual <G> - G and 

' ' meas 

values evaluated from the linear equation. 

(3) Sensitivity of y to errors of the linear equation is given by 

Y = (DGAM/DG) ❖ SD 

’'error ' 
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Table 2. Calibration Data for x-wire Probe Array. 


, .. VEL 

1.21030E 

1.07906E 

9.33095E 

7.96535E 

6.71927E 

5.29811E 

3.96047E 

2.58724E 

1.19316E 

1.91827E 

3;28003E 

4.73371E 

6.13657E 

7.42294E 

8*83925E 

1.01225E 

1*15125E 


STAND DEV 
E 0**2 

M 


, WIRE 1 
4.04599E 00 
3*97028E 00 
3.88373E 00 
3.79260E 00 
3.701 lOE 00 
3.57423E 00 
3.43806E 00 
3.25336E 00 
2.96599E 00 
3.12629E 00 
3.35071E 00 
3.52301E 00 
3.65531E 00 
3.75814E 00 
3.85291E 00 
.3.9.3376E 00 
4.01068E 00 


4.37815E-01 
4.74595E 00 
4.49572E-01 


WIRE 2 
4.38013E 00 
4.30451E 00 
4.22596E 00 
4.14354E 00 
4. 05441 E 00 
3.94312E 00 
3.80764E 00 
3.63442E 00 
3.37164E 00 
3.53613E 00 
3.73504E 00 
3.88813E 00 
4.00981E 00 
4.10600E 00 
4.20058E 00 
4.2725.7E ,00 
4.34596E 00 

5.50605E-01 
7.27388E 00 
4.49399E-01 


Note; Stand. Dev (standard deviation) is based upon V ^ from 

E = + K V and from pressure transducer. 

2 

m. and from least squares routine. 
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VEL 

120 • 
110 . 
100 . 
90 . 
B0 « 
70 . 
60 . 
50 . 
40 . 
30 . 
20 . 
10 . 


- 40 . 

- 35 . 

- 30 . 

- 25 . 

- 20 . 

0.3 

-l.O 

0.4 

- 0.4 

1.2 

3.6 

1.2 

1.7 

0.6 

2. 1 

— 0.8 

- 0.4 

1.5 

0.2 

0.5 

0.4 

1.0 

- 0.7 

Q.e 

1.3 

- 0.3 

0.9 

- 0.6 

- 0.4 

0.9 

- 0.4 

- 0.4 

- 0.2 

0.1 

-O.l 

- 0.2 

- 1.3 

- 0.7 

- 1.2 

0.0 

- 1.9 

- 2.0 

- 0.1 

0.2 

0.0 

- 1.2 

- 0.3 

- 0.9 

0.5 

0.1 

1.0 

0.5 

0.0 

- 0.2 

-O.l 

1.5 

1.8 

2.1 

1.1 

2.2 

4.2 

4.9 

3.9 

7.0 

6.4 


- 15 . 

1.1 

1.1 
2.0 
1.6 

1.4 
0.2 
0.6 
0.1 
0.7 
0.3 
3.1 

6.4 




1.0 

1.9 

2.0 

1.6 

1.3 

0.9 

0.1 

0.0 

0.6 

0.2 

2.5 

6.1 


0.6 

0.4 

0.4 

0.1 

1.7 

2.0 

1.3 

1.6 

1.2 

1.4 

0.8 

0.9 

1.5 

1.4 

0.2 

0.2 

1*0 

0.0 

0.7 

0.4 

0.5 

- 0.2 

0.4 

0.1 

0.2 

0.5 

0.0 

- 0.4 

0.1 

0.6 

- 0.3 

- 0.3 

0.5 

0.2 

0.0 

0.0 

1.0 

0.8 

0.7 

- 0.2 

2.2 

2.9 

1.7 

2.0 

7.7 

7.0 

9.6 

7.6 


O.l 

1.0 

1*0 

0.6 

0.0 

- 0.6 

- 0.2 

- 0.6 

- 0.4 

0.1 

0.8 

7.6 


2 }* 

0 *0 
0.8 
0.2 
“ 0.1 
- 0.6 
- 0.6 
- 1.3 
- 1.1 
- 0*8 
-0.9 
0.1 
4*3 


25 . 

- 0.5 

0.4 

- 0.2 

- 0.8 

- 0.9 

- 1.3 

- 1.6 

- 2.0 

- 1.8 

- 1.6 

- 0.9 

3.1 


30 . 

- 0.8 

0.1 

- 0.6 

- 0.9 

- 1.3 

- 1.6 

-1.9 

- 2.1 

-2.7 

- 2.3 

- 1.1 

2.9 


35. 

0.0 

0.8 

0.1 

- 0.6 

-0.7 

- 1.0 

-1.3 

- 1.6 

- 1.9 

- 1.8 

-1.5 

4.4 


40 . 

1.4 

2.2 

1.5 
1.2 
0.7 
0.4 
0.0 

- 0.8 

“ 1.2 

-1.3 

-l.O 

5.2 


VEL 

- 40 . 

- 35 . 

- 30 . 

- 25 . 

• 

o 

CM 

120 . 

- 0.4 

1.3 

- 0.5 

0.5 

- 1.4 

110 * 

- 4.1 

- 1.4 

- 1.9 

- 0.7 

- 2.4 

100 . 

0.8 

0.4 

- 1.5 

- 0.2 

- 0*5 

90 . 

- 0.4 

- 0.9 

0.6 

- 0.7 

- 1.2 

80 . 

0.3 

- 0.7 

0.5 

0.3 

- 0.7 

70 . 

0.2 

0.2 

0.1 

- 0.1 

0.0 

60 . 

0.1 

0.6 

0.4 

0.7 

0.0 

50 . 

0.9 

1.0 

0.0 

- 0.1 

0.0 

40 . 

0.4 

0.1 

0.3 

- 0.2 

0.0 

30 . 

- 0.3 

- 0.1 

0.0 

0.0 

0.0 

20 . 

- 0.3 

- 0.4 

- 0.4 

- 0.2 

- 0.4 

10 . 

- 0.5 

- 0.5 

- 0.5 

- 0.7 

- 0.8 


- 15 . 

DIFFERENCE 
- 10 . - 5 . 

; ON V 

0 . 

< FPS ) 

5 . 

10 . 

- 1.3 

- 1.3 

- 0.7 

- 0.5 

- 0.5 

- 0.2 

- 1.2 

- 2.2 

- 1.9 

- 2.2 

— 1.5 

— 1.0 

- 2*0 

- 2*0 

- 1.2 

- 1.4 

- 0.8 

- 0.9 

- 1.4 

- 1*4 

- 1.3 

- 1.2 

- 0.1 

— 0 . 2 

- l.l 

- 1.0 

- 0.8 

0.0 

- 0.6 

- 0.3 

- 0.2 

- 0.6 

- 0.3 

0.1 

- 0.3 

0.0 

- 0.4 

0.0 

- 0.1 

- 0.3 

0.0 

0.2 

0.0 

0.0 

0.0 

- 0.3 

0.1 

0.1 

- 0.2 

- 0.2 

- 0.2 

0.0 

0.0 

0.0 

- 0.1 

0.0 

- 0.3 

- 0.2 

- 0.2 

0.0 

- 0.6 

- 0.5 

- 0.4 

- 0.5 

- 0.3 


- 0.8 

- 0.7 

- 0.9 

- 0.8 

- 1.0 

- 0.8 


15 . 

20 . 

25 . 

30 . 

35 . 

40 . 

- 0.2 

0.0 

0*6 

1*0 

0.0 

- 1.7 

- 1.2 

- 0.9 

- 0.4 

- 0.1 

- 0.9 

— 2.5 

- 1.0 

- 0.2 

0.2 

0*6 

- 0.1 

- 1.5 

- 0.6 

0 * 1 

0.7 

0*8 

' 0.5 

- 1.1 

0.0 

0.5 

0*7 

1.0 

0.6 

- 0.5 

0.4 

0.4 

0*9 

1.2 

0.7 

- 0.2 

O.l 

0.7 

0.9 

1.1 

0.7 

0.0 

0.3 

0.5 

1.0 

1*1 

0.8 

0.4 

0.1 

0.3 

0.7 

1.1 

0.7 

0.4 

0.0 

0.3 

0.5 

0.7 

0.5 

0.4 

- 0.1 

0.0 

0.1 

0.2 

0.3 

0.2 

- 0.8 

- 0.5 

- 0.3 

- 0.3 

- 0.5 

- 0.5 


Ui 

O 



VEL 

- 40 . 

- 35 . 

- 30 . 

- 25 . 

— 20 . 

- 15 . 

PERCENT ERROR ON 
- 10 . - 5 . 0 . 

GAMMA 

5 * 

10 . 

15 . 

20 . 

25 . 

30 . 

35 . 

40 . 

120 . 

2.3 

- 4.3 

- 4.9 

- 3.3 

2.3 

7.6 

11.8 

11.8 

11.8 

11.8 

7.6 

2.3 

- 0.7 

- 4.9 

- 4.3 

0.7 

110 . 

2.3 

- 4.3 

- 4.9 

- 3.3 

2.3 

3.4 

ll.S 

11.8 

11*8 




- 3.3 

- 4.9 

- 4.3 

0 .7 

100 . 

0.7 

- 4.3 

- 4.9 

- 3.3 

- 0.7 

7.6 

11.8 

11*8 

11.8 

5.5 

3.4 

2.3 

- 3.3 


- 4.3 

0.7 

90 . 

0.7 

- 4.3 

- 7.0 

- 3.3 

2.3 

7.6 

11.8 

24.4 

11.8 

11.8 

3.4 

2.3 

- 0.7 

— 2 . 8 

—2 * 5 

0*7 

80 . 

0.7 

- 4.3 

- 7.0 

- 5.8 

2.3 

7.6 

11.8 

11.8 

- 0.7 

5.5 


2.3 

- 3.3 

- 4.9 

- 4.3 

0.7 

70 . 

0.7 

- 4.3 

- 7.0 

- 3.3 

- 0.7 

3.4 

ll.B 

11*8 

11.8 

§•5 


2 . 3 

- 0.7 

- 4.9 

- 4.3 

0.7 

60 . 

0.7 

- 6.1 

- 7.0 

- 5.8 

- 0.7 

7.6 

5.5 

n.8 

11.8 

- 0.7 

5.5 



- 3.3 

- 4.9 

- 4.3 

0.7 

50 . 

0.7 

- 6.1 

- 7.0 

- 3.3 

- 0.7 

3.4 

5.5 

- 0.7 

5.5 


- 0.7 

- 0*7 

- 4.9 

- 4.3 

2.3 

40 . 

0.7 

- 6.1 

- 9.1 

- 3.3 

- 0.7 

3.4 

11.8 

11.8 

11.8 

§•5 

3.4 

- 0.7 

- 3.3 

- 0.7 

— 2 . 5 

3.9 

30 . 

0.7 

- 4.3 

- 7.0 

- 5.8 

- 3*9 

3.4 

5.5 

18.8 

- 0.7 

5.6 


- 0.7 

- 3.3 

- 2.8 

- 2.5 

5 .5 

20 . 

2.3 

- 2.5 

- 7.0 

- 8.3 

5.5 

3.4 

5.5 

- 0.7 

11.8 

- 0.7 

3.4 

2.3 

- 0.7 

1.3 

2.6 

10.2 

10 . 

5.5 

- 2.5 

- 7.0 

- 5.8 

- 7.0 

- 0.7 

5.5 

49.6 

- 0.7 

5.5 

3.4 

- 0.7 

- 0.7 

— 2 . 8 

4.6 

13.3 


VEL 

- 40 . 

- 35 . 

- 30 . 

- 25 . 

- 20 . 

- 15 . 

DIFFERENCE 
- 10 . - 5 . 

ON GAMMA < DEG ) 

0 . 5 . 10 . 

15 . 

20 . 

25 . 

30 * 

35 . 

40 . 

120 . 

0.9 

- 1.5 

- 1.4 

- 0.8 

0.4 

1.1 

l.l 

0*5 

0.0 

- 0.5 

- l.l 

- 1.1 

- 0*4 

0.1 

1.4 

1.5 


no . 

0.9 

- 1.5 

- 1.4 

- 0.8 

0.4 

0.5 

l.l 

0.5 

0.6 

- 0.5 

- 0.5 

- l.l 

- 0.4 

0*8 

1.4 

1.5 

- 0.3 

100 . 

0.3 

- 1.5 

- 1.4 

- 0.8 

- 0.1 

1.1 

1.1 

0.5 

0.6 

- 0.5 

- 0.5 

- 0.5 

- 0.4 

0,8 

1.4 

1.5 

- 0.3 

90 . 

0.3 

- 1.5 

- 2.1 

- 0.8 

0.4 

1.1 

1.1 

1.2 

0.6 

- 0.5 

- 1.1 

- 0.5 

- 0.4 

0.1 

0.8 

0.9 

- 0.3 

80 . 

0.3 

- 1.5 

- 2.1 

- 1.4 

0.4 

1.1 

1.1 

0.5 

0.0 

0.0 

- 0.5 

- 1.1 

- 0.4 

0.0 

1.4 

1*5 

- 0.3 

70 . 

0.3 

- 1.5 

- 2.1 

- 0.8 

- 0.1 

0.5 

1.1 

0.5 

0.0 

- 0.5 

- 0.5 

- 1.1 

- 0.4 

0.1 

1.4 

1.5 

- 0.3 

60 . 

0.3 

- 2.1 

- 2.1 

- 1.4 

- 0.1 

1.1 

0.5 

0.5 

0.6 

0*0 

- 0.5 

- 0.5 

- 0.4 

0.8 

1.4 

1.5 

- 0.3 

50 . 

0.3 

- 2.1 

- 2.1 

- 0.8 

- 0.1 

0.5 

0.5 

0.5 

0.6 

0.0 

- 0.5 

- 1.1 

O.l 

0.1 

1.4 

1.5 

-0.9 

40 . 

0.3 

- 2.1 

- 2.7 

- 0.8 

- 0.1 

0.5 

1.1 

0.5 

0.6 

- 0.5 

- 0.5 

- 0.5 

0.1 

0.8 

0.2 

0.9 

-1 .5 

30 . 

0.3 

- 1.5 

- 2.1 

- 1.4 

- 0.7 

0.5 

0.5 

0.5 

0.6 

0.0 

- 0.5 

- 0.5 

0.1 

0.8 

0.8 

0.9 

— 2.2 

20 . 

0.9 

- 0.9 

- 2.1 

- 2.0 

1.1 

0.5 

0.5 

0.0 

0.6 

- 0.5 

0.0 

- 0.5 

- 0,4 

0.1 

- 0.3 

-0.9 

“ 4.0 

10 . 

• 2.2 

- 0.9 

- 2.1 

- 1.4 

- 1.4 

- 0.1 

0.5 

2.4 

1.2 

0.0 

- 0.5 

- 0.5 

0.1 

O.l 

0.8 

- 1.6 

- 5.3 


. Ability of interpolation formulae to reproduce calibration data. ^ 

Note: Data columns represent y values over the range: - 40^ ^ Y ^ 40 . 


Table 3 



